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We present a search for the decay of a B'^ or B'^ meson to a K*^K^ or K*^K''^ final state, using a 
sample of approximately 232 million BB events collected with the BABMi detector at the PEP-II 
asymmetric energy e"'"e~ collider at SLAC. The measured branching fraction is B{B° — > K*°K°) + 
K"^K°) = (0.2l!; j tl\) X 10"^. We obtain the following upper limit for the branching 
fraction at 90% confidence level: 5(5" ^ K*'^ + B{B° K'-"K°) < 1.9 x 10"^ Wc use 
our result to constrain the Standard Model prediction for the deviation of the CP asymmetry in 
B° -i- (j)K° from sin2/3. 



PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 

I. INTRODUCTION 

This paper describes a search for the decay of a 5° 
or B° meson to a K*^K° or K*°K° final state. Hence- 
forth, we use B° K*°K° to refer to both B° and B° 
decays and to the K*'^K° and K*^K'^ decay channels. 
In the Standard Model (SM), K*^K° decays are 

described by 6 — + dss diagrams such as those shown in 
Fig. n Figure ^a) illustrates b ^ d "penguin" transi- 
tions. A so-called rescattering process, effectively a tree- 
level b dvM weak decay followed by the long distance 
production of a ss pair, is shown in Fig. ^b). Other 
rescattering diagrams, e.g., with an intermediate quark 
loop rather than a u quark loop, are also possible. Note 
that the rescattering diagrams can be considered to be 
the long distance components of the corresponding pen- 
guin diagrams, in which the quark in the intermediate 
loop approaches its mass shell. 

The SM prediction for the branching fraction of B^ 
K*^K° is about 0.5 x 10^6 fl-Q. Extensions to the SM 
can yield significantly larger branching fractions, how- 
ever. For example, models incorporating supcrsymme- 
try with R-parity violating interactions predict branching 
fractions as large as about 8 x 10~^ The event rates 
corresponding to this latter prediction are well within 
present experimental sensitivity. Currently, there are no 
experimental results for B'^ — > K*^K^ . Searches for the 
related non-resonant decay _B" K^tt^K^ are reported 
in Ref. H. 

At present, little experimental information is available 
for 6 ^ d transitions. Such processes can provide im- 
portant tests of the quark-flavor sector of the SM as 
discussed, for example, in Ref. j5||. Our study can also 
help to clarify issues concerning potential differences be- 
tween determinations of sin2/3 from tree- and penguin- 
dominated processes, where (3 is an angle of the Unitarity 
Triangle. Such differences can provide a signal for physics 
beyond the SM In particular, our study is relevant 
for the interpretation of the time dependent CP asym- 
metry obtained from B^ — > (j)K'^ decays. (For a review 
of the Unitarity Triangle and sin2/3 measurements based 
on B° (j)K° decays, see Sec. 12 of Ref. 0.) In the SM, 
this decay is dominated by the 6 — > s penguin diagrams 
shown in Fig. |21[a). In addition, sub-dominant SM pro- 
cesses with a different weak phase, such as those shown 
in Figs. |2Jb) and (c) involving the CKM matrix element 
Vub, contribute at a level that is believed to be small 
The deviation of the CP asymmetry in _B" (f)K^ de- 
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FIG. 1: Feynman diagrams for 5° K*'^K'^: (a) penguin 
diagrams and (b) b —> u rescattering diagram. 

cays from sin2/3 because of these sub-dominant processes 
is referred to as AS^o. 

Grossman et al. [9j introduced a method to obtain 
a SM bound on AS^ko, using SU(3) flavor symme- 
try to relate sub-dominant terms such as those shown 
in Figs. 121b) and (c) to the corresponding terms in 
strangeness-conserving processes such as those shown in 
Fig. ^ To determine this bound, measurements of the 
branching fractions of 11 B'^ decay channels are required 
{K*"K°, K*°K°, and hh' with h = cj), uj or and h' = r?, 
r]' or tt'^). Experimental results are currently available for 
all these channels except the two in our study: K*^K^ 
and K*'^K^ . Our measurements will therefore enable this 
bound on AS^^o to be determined for the first time. 
Note that there are not statistically significant signals for 
any of the nine channels for which results are currently 
available. 

Our results might also help to constrain predictions 
for other charmless, strangeness-conserving decays such 
as B'^ —> pTT, in which a dd or uu pair couples to the gluon 
in Fig.nfa) rather than a ss pair (see, e.g.. Table III of 
Ref. 11). 

II. THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the 
BABAR detector at the PEP-II asymmetric e+e~ stor- 
age ring. The data sample consists of an integrated lu- 
minosity of 210 fb~^ recorded at the T(4S') resonance 
with a center-of-mass (CM) energy of y/s — 10.58 GeV, 
corresponding to (232 ± 2) x 10^ BB events. A data 
sample of 21.6 fb^^ with a CM energy 40 MeV below 
the T{4:S) resonance is used to study background con- 



7 




d d 



(C) 

FIG. 2: (a)_CKM Dominant and (b,c) CKM suppressed dia- 
grams for B° (l>K°. 

tributions from light quark e+e" ~^ qq {q = u,d,s or c) 
continuum events. 

The BABAR detector is described in detail else- 
where Charged particles are reconstructed using 
a five-layer silicon vertex tracker (SVT) and a 40-layer 
drift chamber (DCH) immersed in a 1.5 T magnetic 
field. Charged pions and kaons are identified (parti- 
cle identification) with likelihoods for particle hypothe- 
ses constructed from specific energy loss measurements 
in the SVT and DCH and from Cherenkov radiation 
angles measured in the detector of internally refiected 
Cherenkov light. Photons are reconstructed in the elec- 
tromagnetic calorimeter. Muon and neutral hadron iden- 
tification are performed with the instrumented flux re- 
turn. 

Monte Carlo (MC) events are used to determine signal 
and background characteristics, optimize selection crite- 
ria, and evaluate efficiencies. B'^B^ and B^B~ events, 
and continuum events, are simulated with the EvtGen [TH 
and Jetset event generators, respectively. The effec- 
tive integrated luminosity of the MC samples is at least 
four times larger than that of the data for the B^B^ and 
B~^B~ samples, and about 1.5 times that of the data for 
the continuum samples. In addition, separate samples 
of specific B^S" decay channels are studied for the pur- 
poses of background evaluation (see, e.g., the channels 
mentioned in Sec. IIII B|) . All MC samples include simu- 



lation of the BABAR detector response 



III. ANALYSIS METHOD 
A. EVENT SELECTION 

B" K*'^K^ event candidates are selected through 
identification of K*° K+n- and tt+t:- 
decays. Throughout this paper, the charge conjugate 
channels are implied unless otherwise noted. As the first 
step in the selection process, we identify events with at 
least five charged tracks and less than 20 GeV of total 
energy. candidates are formed by combining all op- 
positely charged pairs of tracks, by fitting the two tracks 
to a common vertex, and by requiring the pair to have a 
fitted invariant mass within 0.025 GeV/c^ of the nominal 
Kg mass assuming the two particles to be pions. The 
candidate is combined in a vertex fit with two other op- 
positely charged tracks, associated with the K*''^ decay, 
to form a B" candidate. These latter two tracks are each 
required to have a distance of closest approach to the 
e"'"e~ collision point of less than 1.5 cm in the plane per- 
pendicular to the beam axis and 10 cm along the beam 
axis. Of the two tracks associated with the decay, 
one is required to be identified as a kaon and the other 
as a pion using the particle identification. Charged kaons 
are identified with an efficiency and purity of about 80% 
and 90%, respectively, averaged over momentum. The 
corresponding values for charged pions are 90% and 80%. 

Our study utilizes an extended maximum likelihood 
(ML) technique to determine the number of signal and 
background events fSec. IIII^ . The fitted experimental 
variables are AE, rriEs, and the mass of the K*'^ can- 
di date Mk+^~ , with AE = - Keain ^nd mEs = 
-J^beam - P*b' .^5' where E*g and are the CM en- 
ergy and momentum of the B^ candidate and B^^^^j is 
half the CM energy. M^+tt- is determined by fitting the 
tracks from the K*'^ candidate to a common vertex. We 
require events entering the ML fit to satisfy the following 
restrictions: 

• |AB| < 0.15 GeV, 

• 5.2 < TOES < 5.3GeV/c2, 

• 0.72 < Mk+^- < 1.20GeV/c2. 

Note that virtually all well reconstructed signal events 
satisfy these criteria. 

We further impose the following criteria. The selection 
values are optimized to minimize the estimated upper 
limit on the B" — + K*'^K^ branching fraction by com- 
paring the number of expected signal 2] and background 
events as the selection values are changed. 

• B° criteria: The probability of the fitted B'' 
vertex is required to exceed 0.003. 
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• K*" criteria: K*" candidates are required to satisfy 
I cos I > 0.50, where 0h is the hehcity angle in 
the K*'^ rest frame, defined as the angle between 
the direction of the boost from the i?° rest frame 
and the momentum. 

• K'^ criteria: The probability of the fitted K'^ ver- 
tex is required to exceed 0.06. The fitted K'^ mass 
is required to lie within 10.5 MeV/c^ of the peak of 
the reconstructed K'^ mass distribution. (For pur- 
poses of comparison, one standard deviation of the 
K° mass resolution is about 3 MeV/c^) The K*^ 
decay length significance, defined by the distance 
between the K*° and K'^ decay vertices divided by 
the uncertainty on that quantity, is required to be 
larger than 3. The angle between the K'^ flight di- 
rection and its momentum vector, O^o, is required 
to satisfy cos^^o > 0.997, where the K'^ flight di- 
rection is defined by the direction between the K*'^ 
and K'^ decay vertices. 

• Event shape criteria: To separate signal events 
from the continuum background, we apply selec- 
tion requirements on global momentum properties. 

mesons in T(4S') decays are produced almost 
at rest. Therefore, the i?° decay products are es- 
sentially isotropic in the event CM. In contrast, 
continuum e+e^ qq events at the T(4S') energy 
are characterized by back-to-back two-jet-like event 
structures because of the relatively small masses of 
hadrons containing u, d, s and quarks. As a means 
to separate signal from continuum background 
events, we calculate the Legendre polynomial-like 
terms Lq and L2 defined by Lq = o c Pi ^^'^ 
-^2 = Sr.o.e. ^: (3cos^ 6^ — 1), whcrc Pi is the mag- 
nitude of the 3-momcntum of a particle and 9i is 
its polar angle with respect to the thrust axis, 
with the latter determined using the candidate 
decay products only. These sums are performed 
over all particles in the event not associated with 
the decay ( "rest-of-event" or r.o.e.). Lq and 
L2 are evaluated in the CM frame. We require 
0.374 Lo - 1.179 L2 > 0.15. The coefficients of Lq 
and L2 are determined with the Fisher discrimi- 
nant method . To further reduce the continuum 
background, we also require | cos6'x| < 0.55, where 
6t is the angle between the momentum of the _B" 
candidate and the thrust axis, evaluated in the CM 
frame, with the thrust axis in this case determined 
using all particles in the event except those associ- 
ated with the B^ candidate. 

After applying the above criteria, 3.8% of the selected 
events are found to contain more than one B^ candidate. 
For these events, only the candidate with the largest B^ 
vertex fit probability is retained. 

Our selection procedure eliminates 99.78% and 99.97% 
of the BB and continuum background MC events, re- 
spectively, while retaining 9.8 ±0.1% of the signal MC 



events. 



B. BACKGROUND EVALUATION 

To identify residual backgrounds from B decays, we 
examine and B+B' MC events that satisfy the 

selection criteria of Sec. 1111 Al and that fall within the 
expected signal region of the rngs distribution, defined by 
5.271 < mEs < 5.286 GeV/c^. The events so-identified 
are divided into four categories. 

1. Events containing B^ decays with the same 
K-KTrTT final state as the signal, such as B^ — > 
D'^K^iD^ n^K"^), B° -> D^^i (£)T ^ 
i4:±7r=F7r=F), or B° -> K^tt'^K'^. These channels 
are expected to peak in the signal regions of toes 
and AE but not in the signal region of Mjf+^-. 
The largest number of background events in this 
category arises from B'^ D^K^ (D^ tt^K'^). 
To reduce the contributions of this channel, we ap- 
ply a veto on the n^K'^ mass M^^o based on the 
invariant mass of the and the pion used to re- 
construct the K*°. A veto with 1.813 < M^^o < 
1.925 GeV/c^ (corresponding to ±7 standard devi- 
ations of a Gaussian fit to the M^^o MC distribu- 
tion) removes 64 ± 1% of the D'^K^ background 
MC events but only 4.4 ± 0.6% of the signal MC 
events, where the uncertainties are statistical. Note 
that the reconstructed M^j^o distribution has non- 
Gaussian tails. 

2. Events containing B'^ decays with a kaon misiden- 
tified as a pion, such as B^ ~* K'^K") 
or B° fK° if K+K-). This category of 
background is expected to peak in the mES signal 
region, but not in the M[(+^- signal region, and to 
exhibit a peak in AE that is negatively displaced 
with respect to the signal peak centered at zero. 
The largest number of events in this category arises 
from B^ — > (i>K2: {<t> — > K'^K^). We apply a veto on 
the K^K^ mass Mx+k- assuming the pion candi- 
date used to reconstruct the if *° to be a kaon. The 
veto requires 1.0098 < Mk+k- < 1.0280 GeV/c^ 
(corresponding to ±2.5 standard deviations of a 
Gaussian fit to the Mk+k- MC distribution). This 
selection requirement eliminates 87 ± 1% of the 
(fyK^s background MC events but only 1.2 ±0.3% 
of the signal MC events. 

3. Events containing B^ decays with a pion misidenti- 
fied as a kaon, such as B° D^ir'f (D± tt^K^) 
or B° p"A'° (p° 7r±7r=F). This category of 
background peaks in the mES signal region but not 
in the M^+^r- signal region and exhibits a peak in 
AE that is positively displaced from zero. 

4. AU remaining and B+B~ MC events that 
do not fall into the three categories listed above. 
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such as B° K*°j {K*° K^tt^^), B° ^ 
D^K^ (D=F ^ tJ.-v^,K°,), or B° ^ i K° (r/ ^ 
p^^). These events are characterized both by parti- 
cle misidentification and an exchange of tracks be- 
tween the B and B decays. This class of events 
does not peak in /S.E. 

Based on scaling to the experimental luminosity, 
1.0 event (rounded to the nearest integer) is expected 
for each of the first three categories, and 54 events for 
the fourth category. 

We also consider potential background from the fol- 
lowing source. 

5. Events with the same Kmrn final state as our sig- 
nal but with a K^tt^ S-wave decay amplitude, ei- 
ther non- resonant or produced, e.g., through B^ — > 
Xo*"(1430)A"0 {K^"{U30) X±7r=F) decays. 
These channels are expected to peak in the signal 
regions of toes and AE but not in the signal region 
of Mk+,-. 

There are no experimental results for 
KS°{USO)Kl Studies _16] of B+ K+tt+tt- found 
a substantial i?+ iirQ°(1430)7r+ resonant component, 
however. To evaluate this potential source of back- 
ground, we generate 5° ^ i4:*0(1430)X° {K*^{UZQ) 
K^TT~) MC events. After applying the criteria de- 
scribed in Sec. IIII Al only 1.4 ± 0.1% of these events re- 
main. More importantly, the interference between the 
if*" (890) and S-wave Kn amplitudes is expected to can- 
cel if the detection efficiency is symmetric in the can- 
didate K*^ cos6'h distribution. Through MC study, we 
verify that our efficiency is symmetric in cos 0h to better 
than about 10%. This allows us to treat potential S-wave 
K^t:^ background as an independent component in the 
ML fit. 



C. FIT PROCEDURE 

An unbinned extended maximum likelihood fit is used 
to determine the number of signal and background events 
in the data. The extended likelihood function C is defined 

by 



N 



£ = exp - ^ rQ 



i=i li=l 



(1) 



where N is the number of observed events and are 
the yields of the seven event categories: signal, contin- 
uum background, and the five BB background categories 
from Sec. IIII B1 The correlations between the three fitted 
observables are found to be small ( < 10% in both signal 
MC and background). Therefore, we define the functions 
Vi to be products of three independent probability den- 
sity functions (PDFs), one for each of AE', TOes, and 
Mk+tt- ■ We account for effects related to residual corre- 
lations between the variables through the bias correction 



and evaluation of systematic uncertainties discussed in 
Sees. HVI and M 

The signal PDFs are defined by a double Gaussian dis- 
tribution for AE, a Crystal Ball function |0| for toes, 
and a Brcit-Wigner function for A/^+^-. The param- 
eters are fixed to values found from fitting signal MC 
events. We verify that the signal MC predictions for 
the AE and toes distributions agree with the measured 
results from B° 4>K^ decays 18] to within the ex- 
perimental statistical uncertainties. The (jjK^^ channel is 
chosen for this purpose because of its similarity to the 
K*^K^^ channel. 

Separate PDFs are determined for the continuum back- 
ground and all five categories of BB background item- 
ized in Sec. IIII Bl The background PDFs are defined by 
combinations of polynomial, Gaussian, ARGUS and 
Breit-Wigner functions fitted to MC events, with the ex- 
ception of the PDFs for the S-wave K^tt^ component for 
which the AE and toes PDFs are set equal to those of the 
signal while the M^+^- PDF is based on the scalar Ktt 
lincshape determined by the LASS Collaboration ^2d\. 
All the fits of PDFs to MC distributions yield values of 
per degree-of-freedom near unity. 

The event yields of the continuum and last two cate- 
gories of BB background from Sec. IIII Bl are allowed to 
vary in the fits, while those of the first three categories 
of BB background are set equal to the expected num- 
bers given in Sec. IIII Bl The PDF shape parameters of 
the continuum events are allowed to vary in the fit, while 
those of the five BB background categories are fixed. 



IV. RESULTS 

We find 682 data events that satisfy the selection crite- 
ria. Application of the ML fit to this sample yields l.Olg g 
signal events and 660 ± 75 continuum events where the 
uncertainties are statistical. These results and those for 
the BB background yields are given in TableQ] Based on 
the SM branching fraction predictions of Ref . ,,2j , 5 sig- 
nal events (rounded to the nearest integer) are expected. 
The number of expected continuum events is 619. The 
statistical uncertainty of the signal yield is defined by the 
change in the number of events required to increase the 
quantity —2 InC by one unit from its minimum value, and 
similarly for the other yields. The statistical significance 
of the result, defined by the square root of the differ- 
ence between the value of — 21n£ for zero signal events 
and at its minimum, expressed in units of the statistical 
uncertainty, is 0.28. 

Figure 121 shows distributions for each of the fitted vari- 
ables. To enhance the visibility of a potential signal, 
events in Fig. |3| are required to satisfy Ci{S)/[£i{S) + 
Ci{B)] > 0.6, where Ci{S) is the likelihood function for 
signal events excluding the PDF of the plotted variable 
i = AE, Toes or Mx+^-, and Ci{B) is the correspond- 
ing term for all background components added together. 
The points with uncertainties show the data. The curves 
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TABLE L Results from the maximum likelihood fit. BB back- 
ground categories 4 and 5 refer to the last two categories of 
background itemized in Sec. IIIIBI The yields for the first 
three BB background categories in Sec. llil51 are fixed to the 
estimated values of 1.0 event each. The uncertainties on the 
yields, fit bias, and efficiencies are statistical. 



Parameter 


Value 


Number of events 


682 


Signal yield 


1 0+**'' 


Continuum background yield 


660 ± 75 


B B background category 4 yield 


IT+lf 

' — 71 


BB background category 5 yield 




ML fit bias (signal bias) 


-0.2 ±0.3 


MC signal efficiency 


9.8 ±0.1% 


(including and </> mass vetos) 




Efficiency corrections 




Kg trackmg 


97.8% 


K*° tracking 


99.0% 


Final-state branching fractions 


23.0% 


Overall detection efficiency 


2.2 ±0.1% 




(0.2«|«i)xl0-« 






Significance with systematics (cr) 


0.26 


90% CL upper limit on 


< 1.9 X lO"*^ 


BiB" K-^K") + B{B° ^ K*''K°) 






S.2 5.21 5.22 5.23 5.24 5.25 5.26 5.27 5.28 5.29 5.3 

m [GeV/c^] 




BaBar 





" 11 










'.u'h. 








m 





1 



0.75 0.8 0.85 0.9 0.95 1 1.05 1.1 1.15 1.2 

M^.^ [GeV/c^] 



FIG. 3: Distributions of A_E, rriES, and M^+.^-. The points 
with uncertainties show the data. The curves show projec- 
tions of the ML fit. A selection requirement on the likelihood 
ratio has been applied as described in the text. The solid 
curve shows the sum of all fitted components, including the 
signal. The dashed curve shows the sum of all background 
components. The dotted curve (barely visible) shows the sig- 
nal component. 



show projections of the ML fit with the likelihood ratio 
restriction imposed. 

We evaluate potential bias in the fitted signal yield by 
applying the ML fit to 250 simulated data samples con- 
structed as described below. The number of continuum 
background events in each sample is derived from a Pois- 
son distribution, with a mean set equal to the number 
of continuum events found in the data, i.e., 660 events. 
We generate Ai?, mEs, and Mi(+^- continuum distribu- 
tions for each sample by randomly sampling the contin- 
uum PDFs using the appropriate number of events for 
each sample. The number of BB background events in 
each sample is determined in the analogous manner for 
each of the five BB background categories separately. 
For the first four categories of BB background (all but 
the scalar Ktt component), the Ai?, tties, and M^+tt- 
distributions are generated by randomly selecting the ap- 
propriate number of events from the corresponding MC 



sample. For the scalar Ki: component, the distributions 
are generated by sampling the PDFs. 

The number of signal events in each simulated sam- 
ple is likewise determined from a Poisson distribution, 
with a mean Nf^^ initially set equal to the fitted sig- 
nal yield Ngig = 1.0. The signal Ai?, mEs, and Mji+^- 
distributions are generated by randomly selecting the ap- 
propriate number of signal MC events for each sample, 
is then adjusted until the mean signal yield from 



The ML fit bias is de- 
and is determined to be 



the 250 samples equals Nsig- 
fined by Nuas = N,,g - N^^ 
—0.2 ±0.3 (stat.) events. Therefore, the corrected signal 
yield is N^ig - Nuas = 1-2 events. 

In our study, we can distinguish K*^K^ from K*'^K° 
events from the sign of the electric charge of the K^. 
However, we do not know the flavor of the B meson (i3° 
or i?°) at decay. Therefore, the observed signal yield is 
related to the sum of the B° K*^K" and 5° ^ K*°K° 
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branching fractions through 

(2) 

where e is the overaU detection efficiency, given by the 
product of the MC signal efficiency and three efficiency 
corrections (Tabled). The K'^ and K*'^ tracking correc- 
tions account for discrepancies between the data and MC 
simulation, while the correction for final-state branching 
fractions accounts for the K'^ if", K'^ tt+tt^, and 
j^*o isT+TT" branching fractions, which are not in- 
corporated into the simulated signal event sample. The 
overall efficiency is e = 2.2%. The factor Ngg in Eq. Q 
is the number of BB events in the initial data sample 
of 210 fb~^. We assume equal decay rates of the T(4S') 
to B°B° and B+B- . 

We find the sum of the branching fractions to be 
B{B° K*^K°)+B{B° ^ K*°K°) = {0.2toi1:o.l) x 
10^^, where the first uncertainty is statistical and the 
second is systematic. The systematic uncertainty is dis- 
cussed in Sec. We determine a Bayesian 90% confi- 
dence level (CL) upper limit assuming a uniform prior 
probability distribution. First, the likelihood function is 
modified to incorporate systematic uncertainties through 
convolution with a Gaussian distribution whose standard 
deviation is set equal to the total systematic uncertainty. 
The 90% CL upper limit is then defined to be the value 
of the branching fraction below which lies 90% of the 
total of the integral of the modified likelihood function 
in the positive branching fraction region. We obtain 
B{B^ K*°K°) + B{B" ^ K*°K°) < 1.9 x 10"^. We 
also use the modified likelihood function to determine 
the significance of our branching fraction result includ- 
ing systematics. This result is listed in Tabled 



V. SYSTEMATIC UNCERTAINTIES 

To evaluate systematic uncertainties, we consider ef- 
fects associated with the ML fit, the BB background 
estimates, the efficiency corrections, the total number of 
BB events, and the K'^ tt+tt^ branching fraction. Ta- 
ble ^provides a summary. 

To estimate the systematic uncertainty related to the 
signal PDFs, we independently vary the 11 parameters 
used to characterize the signal AE, m-Es, and AIj^+t^- 
PDFs. The mean and standard deviation of the cen- 
tral AE Gaussian distribution, and the mean of the toes 
Crystal Ball function, are varied by the statistical uncer- 
tainties found by fitting the corresponding quantities to 
data in a recent study of B^ (f>K^ decays . We vary 
the standard deviation of the mES Crystal Ball function 
to account for observed variations between different run 
periods. The width of the M^+tt- Breit-Wigner function 
is varied by ±0.01 GeV/c^ . The remaining six signal PDF 
parameters are varied by one standard deviation of their 
statistical uncertainties found in the fits to the MC dis- 



tributions (Sec. nil Cp . taking into account correlations 
between parameters. For variations of all 11 parameters, 
the percentage change in the signal yield compared to the 
standard fit is taken as that parameter's contribution to 
the overall uncertainty. The total systematic uncertainty 
associated with the signal PDFs is obtained by adding 
these 11 contributions in quadrature. The largest con- 
tributions are from the variations of the AE mean and 
standard deviation (about 0.3 signal events each). 

The systematic uncertainty attributed to the fit bias 
is defined by adding two terms in quadrature. The first 
term is the statistical uncertainty of this bias (Tabled}. 
The second term is defined by changing the method used 
to determine the bias. Specifically, we evaluate this bias 
by generating the AE, tues, and Mj^+„- distributions 
of the fourth BB background category in Sec. IIII Bl us- 
ing the PDFs rather than sampling MC events, for the 
250 simulated data samples: the difference between the 
results of this method and the standard one allows us 
to assess the effect of residual correlations between the 
variables. The fourth category of BB background events 
is chosen because it dominates the BB background. The 
difference between the corrected mean signal yield and 
the standard result defines the second term. 

To estimate an uncertainty associated with the BB 
background, we vary the assumed numbers of events for 
the three BB background categories for which these num- 
bers are fixed, i.e., the first three background categories 
of Sec. IIII Bl Specifically, we independently vary these 
numbers by +2 and —1 events from their standard val- 
ues of unity, and determine the quadrature sum of the 
resulting changes in the signal yield. 

A systematic uncertainty associated with the presumed 
scalar Ktt lineshape is defined by the difference between 
the signal yield found using the LASS lineshape and a 
uniform (i.e., flat) Ktt mass distribution. 

Systematic uncertainties for the K'^ reconstruction ef- 
ficiency, and for the tracking and particle identifica- 
tion efficiencies of the if"*" and 7r~ used to reconstruct 
the ff*", account for known discrepancies between the 
data and MC simulation for these quantities. Similarly, 
the MC simulation overestimates the number of selected 
events compared to data for values of | cos 9t \ less than 
about 0.9. We assign a 5% systematic uncertainty to 
account for this effect. 

The systematic uncertainty associated with the num- 
ber of BB pairs is determined to be 1.1%. The uncer- 
tainty in the K'^, tt+tt^ branching fraction is taken 
from Ref . . 

The total systematic uncertainty is defined by adding 
the above-described items in quadrature. 



VI. SUMMARY AND DISCUSSION 

In this paper, we present the first experimental results 
for the decay B"(B") — > K*^K^. From a sample of about 
232 milhon BB events, we observe 1.0l3^ ^° ^ K*°K^ 



12 



TABLE II: Summary of systematic uncertainties. 



Systematic effect 



Uncertainty 



ML fit procedure (events) 

Signal PDF parameters 0.5 

Fit bias 0.5 

BB background yields 0.1 

Total uncertainty from ML fit (events) 0.7 



Scalar Kn lineshape (events) 



Total systematic uncertainty for B{ x 10® 



+0.0 
-1.4 



Efficiency corrections (%) 

Kg reconstruction 

K**^ tracking 

Particle identification efficiency 
cos St selection requirement 
Number of BB pairs 
B{K'^s ^ 7r±7r=F) 

Total uncertainty from corrections 



1.4% 
2.8% 
0.8% 
5.0% 
1.1% 
0.1% 

6.1% 



+0.1 
-0.3 



Analogous to Eq. (O, the amplitude A' for 

K*^K^ can be expressed as 

A' = V:,V,db' + V:^,Vudb^. (6) 

In contrast to Eq. neither term in Eq. Q is sup- 
pressed by CKM factors relative to the other. As an ef- 
fective tree-level process, it is therefore possible that the 
diagram of Fig. ^b) dominates K*^K'^ decays. 

(This assumption yields the most conservative limit on 

The method of Grossman et al. Q consists of using 
SU(3) flavor symmetry to relate h'^ and 6" in Eq. © to 
a'^ and a" in Eq. ^ to obtain a bound on the quantity 
^^K° defined by 



Vud 



A 



(7) 



with A given by Eq. Q. The bound on ^^i^o is de- 
rived using the branching fractions of 11 strangeness- 
conserving charmless decays: 



< 



0.5- 



2 [B{K*°K^) + B{K*^K^) 



l_B{f^)_ 

B{(l)KO) 



(8) 



event candidates. The corresponding measured sum 
of branching fractions is B{B° K*°K°) + B{B° 
K*^K") = {0.2t°ol t°o:l) X 10-6. We obtain a 90% con- 
fiden£e level upper hmit of B{B° K*°K°) + B{B° 
K*"K°) < 1.9 X 10"*^. This result constrains certain 
extensions of the SM, such as the R-parity violating su- 
persymmetry models described in Ref. 

Our result also can be used to determine an upper 
bound on AS'^^o , as mentioned in the introduction. The 
amplitude A for B^ — <j)K^ can be expressed as 



A 



'* and a" 



(3) 



with — p'^ — p and a" = p^ — p , where p^ is the 
hadronic amplitude of the penguin diagram with interme- 
diate quark i ~ u, qit [see Figs.|3a) and (b)]. The CKM 
factor multiplying a" in Eq. Q is suppressed by ^^(A^ 
relative to the factor multiplying a^, where A = 0.224 
is the sine of the Cabibbo angle. Therefore, the diagrams 
in Fig.|2fa) are expected to dominate i?° 4)K^ decays. 
As described in Ref. 0], AS^^" is given by 



AS^x° — 2 cos 2/3 sin 7 cos 5 \£,4,k'> 



with 



(4) 



(5) 



where 5 and 7 are the strong and weak phase differences, 
respectively, between a" and a'^. 



where the Ci are SU(3) coefficients and where the nine 
final states fi — hh' are specified in the introduction. 
I^if is related to f ^a'" through |^ 



\S.ct>K» I 



+ \^^Ko\^ + 2cos7i?e (iTT^C^ifo 



l + k^Kol +2cos7i?e(^, 



(9) 

The observed rates of strangeness-conserving processes, 
potentially dominated by 6 — > u rescattering transitions 
such as are illustrated in Fig. Q^b), are therefore used 
to set limits on the contributions of the SM-suppressed 
b ^ u terms shown in Figs. |2tb) and (c), i.e., to set 
limits on transitions which cause a deviation of the CP 
asymmetry in i3° — > (l)K^ decays from sin2/3. 

We evaluate a 90% CL upper limit on \AS^jio \ by gen- 
erating hypothetical sets of branching fractions for the 11 
required SU(3)-related decays: K*°K°, K*°K^, and hh' . 
Branching fraction values are chosen using bifurcated 
Gaussian probability distribution functions with means 
and bifurcated widths set equal to the measured branch- 
ing fractions and asymmetric uncertainties. For the mea- 
surements of the branching fractions of the nine channels 
not included in the present study, see Refs. 22, 23]. Neg- 
ative generated branching fractions are discarded. For 
each set of hypothetical branching fractions, we com- 
pute a bound on \ASfpK°\ using Eqs. Q and For 
the unknown phase term cos i5 in Eq. (0} , we sample a 
uniform distribution between —1 and 1. Similarly, the 
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weak phase angle 7 is chosen by selecting values from a 
uniform distribution between 38 and 79 degrees, corre- 
sponding to the 95% confidence level interval for 7 given 
in Rcf. 24|. (A flat distribution is chosen for 7 because 
the likelihood curve in Ref. is non-Gaussian.) For (3, 
we use sin2/3 — 0.687 jl^. For each iteration of variables, 
Eq. lO is solved numerically for I^^Koj. 

We find that 90% of the hypothetical lAS'^^^oj bounds 
lie below 0.42. Our study thus allows the SU(3) bound 
from Ref. 0, viz., lAS'^^o] < 0.42 at 90% CL, to be 
determined for the first time. To assess the contribu- 
tion of the K*^K'^ channel on this result, we repeat the 
procedure described in the previous paragraph with the 
K*'^K'^ branching fraction and uncertainties set to 
zero: the corresponding result is 0.32. Potential future 
measurements of K*^K^ yielding a significantly 

smaller UL and uncertainties would therefore have a sig- 
nificant impact on the jAS'^^oj bound. As a cross check, 
we also determine the SU(3) bound assuming the weak 
phase angle 7 to be distributed according to a Gaussian 
distribution with a mean of 58.5° and a standard devia- 
tion of 5.8° H: this yields |AS'^jfo| < 0.43 at 90% CL. 
Our analysis does not account for SU(3) fiavor breaking 
effects, generally expected to be on the order of 30%. 
However, the method is conservative in that it assumes 
all hadronic amplitudes interfere constructively. 
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